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Abstract 



O , Optical and electrical characteristics of AlGaAs lasers with sep- 

arate confinement heterostructures are modeled by using Synopsys's 
Sentaurus TCAD, and open source software. We propose a modified 
►-j^ , exponential I —V dependence to describe electrical properties. A sim- 

00 I pie analytical, phenomenological model is found to describe optical 

^^ . efficiency, rj, with a high accuracy, by using two parameters only. A 

—1. . link is shown between differential electrical resistivity r = dU/dl just 

1^ I above the lasing offset voltage, and the functional ry(C/) dependence. 

o 



1 Introduction 



H I Alferov |llj, et al., proposed creating semiconductor-based lasers comprising 

the use of a geometrically-narrow active recombination region where photon 
generation occurs, with waveguides around improving the gain to loss ratio 
(separate confinement heterostructures; SCH). That idea dominated largely 
the field of optoelectronics development in the past years. Due to the relative 
simplicity and perfection of technology, solid solutions of AlxGai-^^s are 
commonly used as wide-gap semiconductors in SCH lasers. 



Reaching the threshold current density of these lasers less than IkA/cm'^ 
at room temperature has opened up prospects for their practical application 
and served as a turning point in their production. Now, they are mostly 
used for pumping solid state Nd : YAG lasers, either for high-power metal- 
lurgical processes or, already, in first field experiments as a highly directional 
source of energy in weapons interceptors. Further progress in that direction 
is associated with optimizing the design of laser diodes and, in particular, in 
improving their optical efficiency as well finding methods of removing excess 
heating released. 

In our earlier works we first were able to find agreement between our cal- 
culations of quantum well energy states and the lasing wavelength observed 
experimentally |j2j. Next [3j, we have shown how to considerably improve 
their electrical and optical parameters by finding the most optimal QW width 
and waveguides widths, and type and level of doping j3] . We compared com- 
puted properties with these of lasers produced by Polyus research institute 
in Moscow [5], [6]. By changing the waveguide profile through introducing a 
gradual change of Al concentration, as well variable doping profiles, we were 
able to decrease significantly the lasing threshold current, increase the slope 
of optical power versus current, and increase optical efficiency. 

We have shown also [7J that the lasing action may not occur at certain 
widths or depths of Quantum Well (QW), and the threshold current as a 
function of these parameters may have discontinuities that occur when the 
most upper quantum well energy values are very close to either conduction 
band or valence band energy offsets. These effects are more pronounced 
at low temperatures, and may be observed also, at certain conditions, in 
temperature dependence of lasing threshold current as well. 

One of the fundamental laser characteristics is their optical efficiency, 
7], the ratio of optical power generated, L, to electrical power supplied, 
7] = L/{U ■ /), as well dependence of rj on current or voltage. We propose 
here a simple analytical, phenomenological model for description of / — l^ 
characteristics near and above lasing offset voltage Uq, and we show that that 
model may be used for description of ri[U) with a high accuracy. At the same 
time we obtain a link between differential electrical resistivity r = dU/dl just 
above Uq, with the functional r){U) dependence, where r is also an important 
experimental characteristic of a laser device. 

For simulations, we use Sentaurus TCAD from Synopsys [8], which is an 
advanced commercial computational environment, a collection of tools for 
performing modeling of electronic devices. 



2 Lasers structure and calibration of modeling. 

We model a laser with lOOOfim cavity length and 100/im width, with doping/ Al- 
content as described in Table [1] 

Synopsys's Sentaurus TCAD is a flexible set of tools used for modeling a 
broad range of technological and physical processes in the world of microelec- 
tronics phenomena. It can be run on Windows and Linux OS. Linux, once 
mastered, offers more ways of an efficient solving of problems by providing a 
large set of open source tools and ergonomic environment for their use, mak- 
ing it our preferred operating system. We find it convenient, for instance, 
to use PerO scripting language for control of batch processing and changing 
parameters of calculations as well for manipulation on text data files, and 
TcO for manipulating (extracting) spacial data from binary TDR files. A 
detailed description, with examples of scripts, is available on our laboratory 
web sitai. 



The results for Ith and S = dL/dl, in this paper, are all shown normalized 
by I^f^ and 5*0, respectively, which are the values of If^ and 5* computed for 
the reference laser described in Table [TJ 

We neglect here the effect of contact resistance, Rx, by not including 
buffer and substrate layers and contacts into calculations (compare with 
structure described in Table [T]). An estimate, based on geometric dimen- 
sions of substrate layers and their microscopic parameters (doping concen- 
tration, carrier mobility) gives us a value of Rx of the order of ISmQ. At 
lasing threshold current of 0.1 A, that small resistance will cause a difference 
between computed by us lasing offset voltage Uq and that measured one by 
about ImV only. We still will however have a noticeable contribution from 
Rx to differential resistance dU/dl. 



3 Methods of data analysis. 

3.1 Threshold current and L{I) dependence. 

The most accurate way of finding Ith is by extrapolating the linear part 
of L{I) to L = just after the current larger than Ith- We used a set of 



"'^Perl stands for Practical Extraction and Report Language: http://www.perl.org 
^ Tool Command Language; http://www.tcl.tk 
^http://www. ostu.ru/units/ltd/zbigniew/synopsys.php 



gnuplot and perl scripts for that that could be run semi-automatically, very 
effectively, on a large collection of datasets. One should only take care that 
the data range for fiting is properly chosen, since L{I) is a linear function in 
a certain range of / values only. The choice of that range may affect accuracy 
of data analysis. 



3.2 Uq from fiting U{I) dependence 

An exponential U{I) dependence is found to work well at voltages which are 
well below the lasing offset voltage Uq- Near the lasing threshold, we observe 
a strong departure from that dependence, and, in particular, for many data 
curves a clear kink in U{I) is observed at Uq. We find that a modified 
exponential dependence describes the data very well: 



/(f/) = Ith ■ expiA .{U-Uo) + B-iU- UoY), for U < U^ ^^^ 
I{U) = Ith ■ exp{C ■{U-Uo) + D-{U- Uof), for U > Uo 

where Ith, Uq, as well A, B, C, and D are certain fiting parameters. Equa- 
tion[T]offers a convenient interpretation of physical meaning of it's parameters 
Ith and Uq: I{Uq) = Ith- 

3.3 Differential resistance 

The above function (Eq. [1]) is continuous at Uq, as it obviously should, but 
it's derivative is usually not. Figure [T] shows a few typical examples of I{U) 
datasets. The lines were computed analytically by using Eq. [H after finding 
all parameters with the least-squares method. 

Since ([T]) may have a discontinuous derivative, using it to find out differ- 
ential resistance at Uq is ambiguous. From Eq. ([1]), at f/ = f/o, we will have 
dU/dl = j-^ on the side U < Uq and dU/dl = j-^ on the side U > Uq. 
Hence, the parameter C may be interpreted in terms of differential resistivity 
just above Uo'- 



^_ (2) 

C-Ith 



We find from data analysis, for instance for the third dataset in Figure 
[H that dU/dl ^ 50mQ, which, together with estimated contact resistance 
Rx = 13mQ gives good qualitative agreement with the differential resistance 
expected for real lasers, where it is in the range of about 50 — 80mQ ([5j and 

El). 



3.4 Doping dependencies 

Figure [2] shows the dependence of parameter D in Eq. [T]on n-, and p-emitters 
doping concentration, for a very broad range of doping concentrations in 
other regions (this is "N-N" type of doping; see description of Table [T]). Due 
to large scatter of the parameters obtained by the least-squares fiting, we 
do not distinguish between datapoints that were obtained for various doping 
concentrations in waveguides or in active region: the dominant factor on 
values oi C OT D parameters is doping concentration in emitter regions. 

We observe also that a correlation exists between values of C and D 
parameters, as illustrated in Figure [31 The line in Figure [2] was obtained by 
using the least-squares fiting method to all the data points displayed there, 
with the following simple function: 



D = -40.073 + 8.6 ■ 10-^ ■ (41.4 - C) 



3.75 



(3) 



It is convenient to rewrite Equation [T] in dimensionless variables. In case 
oi U > Uq we have then: 



i[U 



exp{^-iu-l)-[l + P-Ui-iu-l)]) 



(4) 



where we defined: i{u) = I{U)/Ith and u = U/Uq, a = r ■ Iih/Uq, (3 = 
Uq ■ D/C, and we used also Eq. |2j 

Let us estimate the range of reasonable values of /3 parameter. 

The function [3] would give the ratio D/C -^ +oo for C decreasing to 
(which corresponds to decreasing doping concentration in emitter regions to 
0). That function will pass through at values of C ~ 8.9, which corresponds 
to doping in emitters of around 2 • lO^^cm^^, will have minimum of value 
~ —1.59 at concentrations corresponding to f« 2 • lO^^cm"^, and will increase 
to —1 at larger doping concentrations. The practical range of interest in our 



case is not at the lowest doping concentrations in emitters, since than other 
laser parameters deteriorate. We are left with D /C values that are important 
to us in the range between ?^ — 1.6 and 0. 

Hence, the corresponding range of /3 values that is of our interest is be- 
tween K. —1 and 0. 



4 Optical efficiency 

4.1 A simplified approach 

It is tempting to try a simplified version of HI when the expression under 
exponent is << 1. We have in that case the following approximation on 
optical efficiency: 



„r„l - -^- . (»-l)+/3-(n-l)^ (5) 

'/I"; — Uq-u q:+[(«-1)+/3-(m-1)2] 

Figure H] shows a few example curves computed with Equation [5l The 
accuracy of these results, if compared with real data (not shown on that 
figure), appears to be far from desired; Eq. Omay be treated as a very rough 
approximation only. 



4.2 Exact result 

Let us use however the full version of Equation [T] (for U > Uq), for computing 
optical efficiency. We have then: 



jj(u) = A . ^W-^ (6) 



Uo u-i(u) 



where i{u) is given by Eq. HI 



Figure[S] shows that an excellent agreement is obtained when these analyt- 
ical formulas are used for approximating optical efficiency directly computed 
from modeling data. The parameters used in fiting the data are shown in 
Figure [6l As seen on this figure, the value of /3 for some datapoints is lower 
than —1. This should not be considered contradictory to our estimate of the 
range of possible values of (3: our calculations of r] do not take into account 

6 



the nonlinearity of L{I) dependence, which may be large, especially for low 
doping concentrations in emitters, and that will effectively cause decrease of 
P value. We see also that the parameter a is too small (i.e. {u — l)/a is too 
large) in realistic cases to allow using the simplified Equation [5l 
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Figure 1: Examples of typical I — V characteristics for a few combination 
(as described in the Figure) of doping concentrations (n- and p-emitters 
concentration first, followed by waveguides and active region concentrations). 
We show the narrow region near the lasing threshold, only. The curves are 
computed by using fiting parameters {Ith, A, B, C, D, Uq) of equation [H and 
after that voltage and current are normalized by Uq and Ith-, respectively. 



Table 1: Structure of AlGaAs SCH laser layers used in computer model- 
ing. Values of doping concentrations listed in rows 4-6 of this table are 
approximate only. In this paper we make reference to two kinds of doping 
concentrations in these regions, that we name by "N-N" and "N-P". "N-N" 
means that in two waveguides there is the same N-type doping concentration 
with doping type in active region of n-type as well. "N-P" means" that one 
waveguide is of n-type doping and the other one of p-type doping. Active 
region has then n-type doping. The "real" lasers have N-N type of dop- 
ing, however, setting "N-P" type of doping in computer modeling leads to a 
significant improvement of device characteristics \4\- 



No 


Layer 


Composition 


Doping [cm ^] 


Thickness [fj,m] 


1 


n-substrate 


n-GaAs (100) 


2 ■ IOI8 


350 


2 


n-buffer 


n-GaAs 


1 • 1018 


0.4 


3 


n-emitter 


Alo,5Gao.5As 


1 • 10^8 


1.6 


4 


waveguide 


Alo.ssGao.GjAs 


none (n ^ 10^^) 


0.2 


5 


active region (QW) 


Alo,o8Gao,g2As 


none (n ^ 10^^) 


0.012 


6 


waveguide 


Alo.ssGao.GjAs 


none (n ^ 10^^) 


0.2 


7 


p-emitter 


Alo.5Gao.5As 


1 • 10^8 


1.6 


8 


contact layer 


p-GaAs 


4. 1019 


0.5 




2.5e+18 



p-, n-doping concentration of emitters regions [cm ] 



Figure 2: Dependence of D on n-, and p-emitters doping concentration, for a 
broad range of doping concentrations in other regions: n-type concentration 
in active region changes between 10^^ and lO^^cm"^, while n-type concentra- 
tion in waveguide regions between 10^^ and 10^ 
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Figure 3: Correlation between parameters D and C for a broad range of 
doping concentrations in all regions: n-type concentration in active region 



changes between 10^*^ and 10 
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cm 



-3 



, n-type concentration in waveguide re- 



gions is between 10^^ and lO^^cm "^, and in emitter regions it changes between 



10^^ and 2.5 • lO^^cm^^. 



The line is described in the text. 
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Figure 4: Optical efficiency as a function of voltage, computed by using 
Equation [5l for a few values of parameters a and /3, as shown in the Figure. 
It has been assumed that S = 1.25W/A and Uq = 1.65V^. 
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Figure 5: The lines show optical efficiency computed with equations [6] and HI 
while symbols show optical efficiency directly obtained from modeling data. 
The results are for "N-P" type of doping structure (different type of doping 
for both waveguides), for doping concentration in waveguides of 1 • lO^^cm"^, 
in active region 5 • lO^'^cm"^, and in emitter regions as shown in the Figure. 
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Figure 6: Parameters 5*, Uq, a and /3 as a function of doping concentration 
in emitters, that were used to draw curves in Figure O The hues are to guide 
the eyes, only. 
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5 Summary 

Computer simulations using Sentaurus TCAD from Synopsys were used for 
performing modeling of electrical and optical characteristics of SCH lasers 
based on AlGaAs. 

A modified exponential I — V dependence (Equations [T] and H]) is proposed 
to describe electrical properties. 

That simple analytical, phenomenological model is found to describe one 
of the most fundametal laser characteristics, optical efficiency, f]{U), with 
a high accuracy, by using two parameters only (except of 5* = dL/dl, Ith, 
and Uq). At the same time we obtain a link between differential electrical 
resistivity r = dU/dl just above lasing offset voltage, with the functional 
ri{U) dependence. 

The proposed model is useful for both, analysis of computer modeling 
results as well experimental data on real devices. 
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